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Abstract 

We derive the first £ 2 -correction to the instanton partition functions of J\f = 2 Super Yang-Mills (SYM) in 
four dimensions in the Nekrasov-Shatashvili limit £2 —> 0. In the latter we recall the emergence of the famous 
Thermodynamic Bethe Ansatz-like equation which has been found by Mayer expansion techniques. Here we 
combine efficiently these to field theory arguments. In a nutshell, we find natural and resolutive the introduction 
of a new operator V that distinguishes the singularities within and outside the integration contour of the partition 
function. 
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1 Introduction 


The Omega background was first introduced in N = 2 SUSY gauge theories to regularise the infinite volume of 
M 4 in the computation of instanton contributions to the partition function by localisation [1], More recently, it 
has also proven to be a formidable way to preserve integrability of these theories upon deformation. In fact, this 
background is characterised by two equivariant deformation parameters e\ and £2 associated to the breaking of the 
Lorentz invariant four dimensional space into C x C, but still the Nekrasov instanton partition function exhibits 
an integrable structure in the form of covariance under the Spherical Hecke central (SHc) algebra [2, 3] (which 
is formally equivalent to a algebra). The presence of this algebra shed light on the conjecture by [4, 5] of a 
duality between these four dimensional theories and the family of Toda conformal field theories in two dimensions. 1 
Moreover, the SHc algebra is closely related to a tensorial version of the integrable Calogero-Moser Hamiltonian 
and led to the construction of one of the most basic objects of quantum integrable theories, namely a (instanton) 
/ 1 '-matrix [9], 

In the Nekrasov-Shatashvili (NS) limit £2 — > 0, the theory possesses only one (non-zero) equivarant parameter 
£1 and, besides, its 2D underlying integrable structure is better understood [ 0], The latter defines, indeed, a quan¬ 
tisation of a Hitchin system associated to the Seiberg-Witten curve (see [11] and references therein) which can now 
be characterised by a quantum curve [12]. And this curve is, actually, equivalent to a TQ-relation upon a (quan¬ 
tum) change of variables [13], while the new integrable system is the bi-spectral dual of the previous one [ 4], Of 
course, of great interest for integrable model theory is the TQ-relation, which can also be obtained directly from the 
Nekrasov partition function by extremising the sum over Young diagrams [ 5], a technique that goes back to [16]. 
Moreover, another key feature of integrability has come out of the blue since the original paper [10] in the form 
of a non-linear integral equation (NLIE) in the complex plane [ 1 7] resembling a Thermodynamical Bethe Ansatz 
(TBA) equation [18]. Subsequently, this equation was derived in full detail by [19, 20] upon using the Mayer cluster 
expansion technique [21, 22], 

In the present manuscript, we go beyond this limit and present an explicit formula for the first £ 2 -correction to the 
prepotential. This result should pave the way to a better understanding of the full £2 -deformation and its meaning as 
quantum integrable system characterised by a TBA/NLIE. This should ultimately lead to a richer algebraic structure 
similar to the SHc Hopf algebra as we know how, along the other way around, the NS limit of the SHc algebra can 
be obtained [23]. For convenience, we focus on N = 2 SYM with a single SJJ ( N c ) gauge group and a number of 
fundamental flavours. However, we believe that our mathematical construction and results be easily generalisable to 
arbitrary quivers (along the lines of [23]). 

In fact, with respect to the latter reference (exploiting the methods of Mayer cluster expansion) we shall add an 
efficient field theory argument. The latter is firstly used to disentangle the long- and short-range interactions between 
instantons: this is performed in section 2. Then, we treat the short-range interactions by the Mayer cluster expansion 
technique in section 3 and achieve formulae (3.6) and (3.28) concerning the (short-range) prepotential 2 . Eventually, 
the two arguments are combined in section 4 to build up our proposal for the prepotential at the order 0(e 2 ). In the 
last section, this proposal is checked in the simple case of pure (7(1) SYM up to the order 0(A 4 ). A more detailed 
computation involving only the combinatorics of the Mayer expansion is to appear shortly [24], 

'The conjecture has been partially proved either using the basis of AFLT states [6], or a set of generalized Jack polynomials [7, 8], 

"This part is rather technical and its details could be skipped at first reading. 
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2 Nekrasov partition function as Gaussian correlation of exponential fields 


In the expression derived by N. Nekrasov in [1 ], the gauge coupling expansion of the instanton partition function for 
J\T = 2 SU(N c ) SQCD takes the form of a series over coupled integrals 


Z = 



N 

2—1 


dfa 

2iir 


N 


n *(&-&)> 

i,j =1 
i<j 


(2.1) 


where the ((-background cqinvariant deformation parameters £i, £2 are supposed to have a positive imaginary paid, 
A is the dynamical scale and e + = £1 + £ 2 . In formula (2.1) and then in the following, integrations are performed 
along a contour that surrounds the upper half-plane, including the real axis but excluding possible singularities at 
the complex infinity. The potential Q(x) is a rational function that can be expressed as the ratio of matter and gauge 
polynomials, 


Q(x) 


Ufil (x-rnf) 
Yl^x - cii)(x + e+- ai)' 


( 2 . 2 ) 


where Nf and N c denote respectively the number of flavors (fundamental massive hypermultiplets) and colors 
(adjoint gauge multiplets). The coulomb branch vacuum expectation values (vevs) are assumed to be real but inside 
the integration contour, i.e. a; £ M + iO. The universal kernel for N = 2 SYM with a single gauge group (A \ 
quiver) is 


K[x) 


x 2 (x 2 — £+) 

(x 2 — £ 2 )(x 2 — e 2 ) 


(2.3) 


Since we consider nested integrals, we need a prescription for the poles at d>i = <f>j ± £1 and (fy = cf)j ± £ 2 : we 
assume that the variable cj)j is real, taking the poles at 6 t = 0j + £1 and 6 t = <i> 3 + £ 2 , but not their counterparts at 
0 % = <pj — Ei and 4>i = (j>j — £ 2 . For other gauge groups, the prescription is more subtle and has been discussed in 
[25], 

In order to exploit quantum field theory methods, we factorise the kernel into two pieces 


K(x) = (1 + E 2 p{x))e £2k ( x \ with p(x) 


£2 , / n 2ei T £2 \ 

x ,2 _ e 2 ’ fc ( g ) = £ 2 log ^1 ~ ^2 _ g 2 £ 2 J , 


(2.4) 


where the ’universal’ kernel p(x) contains singularities at x = ±£2 that pinch the integration contour only in the 
NS limit £2 —> 0. It can be interpreted as a strong short-range interaction responsible for the formation of bound 
states of instantons [19, 20]. On the contrary, the long-range interaction k(x) is free of such singularities and can 
be treated by standard Mayer expansion methods [26]. Note that the logarithmic branch cut is not relevant in our 
context since at a finite order in £ 2 , the kernel k{x) is a rational function of x of order 0(1). 

Now, long- and short-range interactions can be separated by the introduction of a Gaussian field X(x) with 
propagator 

(X(x)X(y)) = e 2 k(x - y ). (2.5) 

The corresponding action is a kinetic term with a kernel being the inverse of the propagator, 


5|0ng[X ] “ 2£ 2 / 


dxdy t(x - y)X(x)X(y), f t(x - z)k(z -y)^jz = 2iir5(x - y), 


(2(7r) 5 


2in 


( 2 . 6 ) 


and any correlator 


is understood as Gaussian functional integral average, 

(•■■) = — J— [ DXe~ Sl °"z [x] ■■■ , A/iong = [ DXe~ Slo ^ x] . 
A/long J " J 


(2.7) 
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As a consequence of the Wick theorem or Gaussian integration 3 * , 


N 

e x(< k) 

2=1 

This property allows us to perform an Hubbard-Stratonovich transformation on the instanton partition function (2.1) 

oo n c -N /.tv , , tv 

z = (z short [x]), with z S h 0 Tt[x] = Y n (1 + e 2 P(<l>ij)), (2-ii) 

N=0 ‘ ^ 2=1 Z7r i,3 =1 

i<j 

where we introduced the (finitely) renormalised gauge coupling q = e _ 2 fe ( 0 ) £2 Ae + /ei and the shortcut notation 
4>ij = <t>i — 4>j- Before performing the £ 2 -cxpansion of the quantity Z^ mn X] in the next section, we would like to 
make a remark. It is possible to introduce a second field x( x ) with Gaussian action 

^shortM = * [ ^%r(x-y)x(x)x(y), (2.12) 

Z£2 J [2iir) z 

and the kinetic kernel being the inverse of the logarithm of the propagator, 

/ dz 

t(x - z)k(z - y)^ = 2md(x-y), k(x) = 1 log (1 + e 2 p(x)). (2.13) 

The Nekrasov partition function takes the form of a Gaussian average with respect to the two independent fields, 

Z = ^exp (X/ ' ( 2 - 14 ) 

In this expression the usual Wick regularisation ( e.g . no equal times contractions) will be needed for the logarithmic 
divergence of k(x) in x = 0, as for the Coulomb gas case of 2D conformal field theories. Instead, without pinching, 
i.e. p = 0, then y = 0 and the (combinatorial) results by [26] can be easily obtained by saddle point method (cf. 
also below). This simple fact gave us inspiration for the following methodology. 


_ e ^k(0)£2N 
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( 2 . 10 ) 


3 Treatment of the short range interaction 


In order to derive the expression of Z s ^ on [X] at subleading order in e 2 , we will treat the quantity U(x) = qQ(x)e x(x> 
as a meromorphic potential. This treatment is at the moment not fully justified since X(x) is actually a quantum 
field, but the result we obtain coincide with an alternative, rigorous derivation by performing the Mayer cluster 
expansion of the original partition function Z [24], 

3 Alternatively, this identity can be obtained by simple Gaussian integration of the partition function with a source term 

Ziong[J(x)\ = J DXe~ s '^ [x]+S AAxWd* = 2long [0]e^ S (2.8) 

For this formula is the infinite dimensional d — > oo version of 


(e XlJl e X2j2 


e XdJd ) = v/det T 


n 


dXi -s .T. XiTijXi T. XiX 
—=e ! ' J=1 e i=1 

s/2^ 



(2.9) 


with propagator G — T , where we choose, for the continuum limit of the external field Ji —> J(x), the configuration of point-like sources 
J(x) = 3( x ~~ c f [28] for the connected application to gluon scattering amplitudes ( cf. also the last section. Perspectives). 
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The Mayer expansion [21, 22] of the short-range prepotential /short[A 7- ] = e 2 log Z^faX] is a sum over con¬ 
nected clusters A; with l vertices in which two vertices are connected by at most one link. To each vertex i e VTA;) 
is associated the integration measure Ufa)d())i/2m, and to each link < ij >G /T(A;) the short range interaction 

^2Pfa - <Aj), 


00 ^—(Z—1) 

/"short[A] = E S 7(A0" 
1=1 A, y u 


n n £ 2 p{4>i ~ fa). 

ieV{Ai) <ij>eE(Ai) 


(3-1) 


The symmetry factor <j(Ai) is the cardinal of the group of automorphisms of the cluster A;. It is also useful to 
introduce the dressed vertex K,i noi t(x) obtained as the functional derivative 


lThort(*) = 2mU (x) 


dJ~i short 

5Ufa )' 


This quantity expands as a sum over rooted clusters Af with root x, 


(3.2) 


1 short (-'T 
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(3-3) 


There are two main differences with respect to the prepotential (3.1) both due to the fact that we arc dealing with 
rooted clusters: we do not integrate on the root variable x, and a( Af) < a(A /) as, of course, the automorphisms of 
a rooted cluster keep the root fixed. The short-range free energy /"short [A] can also be expanded over rooted clusters 
due to the property (B.3) of [26], 
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(3-4) 


<ij>£E( Af) 


It was shown in [20] that at first order in £2 the potentials can be transferred to the root of the clusters Af, and the 
remaining integral evaluated exactly to give 4 


1 = 1 J J 

OO 

Y^fa) = J 2 h U W l = ~ M 1 - Ufa)). 
1=1 


The first of these gives rise to the Nekrasov-Shatasvi 1 i action (and thus prepotential as £2 —? 0) once inserted in the 
first of (2.11). Here we want to go beyond this order and give a proposal, which has a little conjectural basis, for the 
sub-leading £2 order. A naive approach would consists in employing a linear approximation for the potential, 

Ufa) ~ Ufa + fa - x)U'(x) + 0(4). (3.7) 


However, this does not work because it overlooks the contributions of the poles of the potential. For instance, the 
contribution of a rooted cluster with a single p-link and two vertices reads 


Ufa) Jpfa — y)u(y)^~ = \u(x)u(x) + = ^ufa) [ufa + E 2 vufa + o(4)\ , 


4 The integral 


has been computed in [20]. 
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where, for any meromorphic function U(x), we define U(x) + = t/ reg . (x + £ 2 ) + [/ s j ng . (x — £ 2 ) or at infinitesimal 
level the operator V 

VU(x) = U' eg ,(x)-U' siDg Xx), (3.9) 

by means of the decomposition U(x) = U re g .(x) + f/ s j ng (x) with U K& _{x) analytic inside the contour of integration 
while C/ging. (x) analytic outside. If U(x) had no singularities inside the integration contour, we could replace the 
symbol V with a partial derivative d x and the result would be the same as in the linear approximation (3.7). In this 
case, a formula for K;hort(x) exact in £2 has been derived in [29] based on an earlier formula obtained by Moore, 
Nekrasov and Shatashvili in [30]. Unfortunately, this case is not relevant to the computation of Z s h ort [X] and we 
present here an alternative method. 

3.1 Articulation links and irreducible clusters 

In the clusters A; it is important to distinguish two types of links. An articulation link of a connected cluster is a link 
whose removal break the cluster into two disconnected pieces [31]. The cluster is said to be (1 particle) reducible. 
Note that if the cluster is a tree, all its links are articulation links. A cluster without any articulation link is said to 
be (1 particle) irreducible. An example of irreducible clusters is the necklace, i.e. a ring of l vertices and l links. In 
[24], we have computed the contributions of all necklaces and observed that they do not exhibit any corrections of 
order 0(e 2 )- We have further computed the irreducible clusters with four vertices for a potential corresponding to 
U( 1) SYM and we have found the same conclusion. It led us to conjecture that irreducible clusters do not provide 
any correction at the subleading order (although higher order corrections are present). 

We need a stronger form of this conjecture. The £ 2 -corrections we are trying to compute are associated to the 
links of the cluster. Heuristically, the total ^-correction is the sum of the corrections brought by each link and 
due to a sort of linearisation. Links of a generic cluster can be either articulation links or not. We claim that only 
articulation links bring non-zero subleading £ 2 -corrections. To justify this claim, we consider a link of an arbitrary 
cluster which is not an articulation link, and we compute its contribution to /short at the order 0(e 2 ). We denote 
the two extremities of the link x and y, and Af’ y the bi-rooted cluster with l vertices obtained by removing this 
link. Since we are interested in the £ 2 -correction of a specific link, it is enough to consider the leading order of the 
integrals associated to the cluster A^’ y . At this order, the potential of each vertex can be transferred to a root, say x, 
and produces U(x) 1 . The remaining integral contribution will be denoted T)(x, y), and the integral associated to the 
original cluster reads 

^ 1 = f -^^ 2 P( x -y) u ( x ) l ^i( x ,y), Fi(x,y) = £j~ l f n ^ II £ 2 p{(j>i-<t>j). (3.10) 

J ^ J iev(A^y) <ij>£E(A*’ y ) 

It is useful to rewrite this expression as 

= J ^ u ( x ) l v i ( x ,x), ui{x,z) = J ^-Fi(x,y)p(y-z), (3.11) 

where the function ui(x, x) would contain possible £ 2 -corrections coming from the p-link we singled out. 

Due to the invariance under translations of theT)(.x, y ) must be a function of the difference (x—y ): T/(.i'. y) = 
T;(x — y), and hence i/;(x, z) = vi(x — z) as well. Moreover, Tj(x) is even as can be shown using the sign flip 
4>i —> Convolution with p-kernels shift the poles in the upper half plane of +£2 and of —£2 in the lower half 

plane. It turns out that T;(x — y) can only have (multiple) poles at x — y = ±n £2 with n = 1, 2,.... In particular it 
has no pole at the complex infinity. We deduce the expansion 
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r l(x - y) = 7 n,m 

n =1 
m =0 
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(x — y ) 2 — n 2 £n 


titT-t 


(3.12) 
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where due to the scale invariance under x, y —> ax , ay, £2 —> ae 2 and d t —>■ ad>i, the coefficients 7 n ,m are C- 
numbers independent of £ 2 . Their value depends on the explicit form of the cluster A'j /y , but at fixed / only a finite 
amount of coefficients are non-zero. Thanks to the formula 
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(3.13) 


it is possible to separate T;(x — y) into regular and singular part with respect to one of its variables, and perform the 
integration over the p -link that has been singled out to obtain i'z(O). We observe that vi{ 0) is a C-number independent 


of £ 2 . 


r'z(O) — 'y ] 7 n,m ( 1) 
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m+ 1 'y ^ 

fc =0 


m + k 
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(2 n(n + 1 )) 
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(3.14) 


This is a strong hint that p -links that are not articulation link do not provide any ^-correction. 


3.2 Considering only trees 


To obtain a better insight on the contributions of articulation links, we first restrict the summation over clusters with 
a tree structure Aj, } and examine the corresponding (tree) dressed vertex 


Y t {x) = U{x)Y j Y, 


_-(*-!) 




1=1 A? 


<7 (A; 


T,U J ieV(AS,.) 
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(3.15) 


<ij>£E( A5,,) 


i^x 


Due to the restriction upon tree clusters, this function obeys the following functional equation [20], 

Y t (x) = U (x) exp (^j p(x - y)Y T (y)^r^J = U(x)e^ T ^ + , 


(3.16) 


the second equality being obtained by an explicit computation of the convolution with the kernel p. Expanding in £2 
with the potential U(x) considered finite, we find 

Y t (x) = Yp\x) + £ 2 Yfi\x) + 0{e |), Y±°\x) = Z 7 (x) e 3 y T 0) (*) Y± l \x) = ^ ^ VTj 0 ) (.r), (3.17) 

2 -Yf, 0) (x) 

with the operator V defined in (3.9). At first order, the generating function can also be expressed in terms of the tree 
function T(z) 5 

Y^\x) = 2T(U(x)/2). (3.19) 

It is enlightening to recover this result from a direct computation. In fact, we have the property that the potential at 
each vertex can be transferred to the root, giving 


Y±°\x) = Y^U(x) 1 Y^ 


1=1 


_-(i-i) 
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d<j>i 


2 m 

iev(A* p <ij>eE(AU 

i^x 


n £ 2 p^-<t>j). 


(3.20) 


s The tree function is related to the Lambert W function on the principal branch as T(z) = — W(—z). It satisfies 
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zmT(z)e- Tl ’\ T{z) = jr ■ 
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Figure 1: A rooted tree decomposed into two rooted subtrees (black and red) by choosing a link (in blue). 


As shown in [20], the remaining integral depends only on the number of links and not on the explicit form of the 
tree. It evaluates to (e 2 /2) 1 ' '. The remaining sum over rooted trees with the appropriate symmetry factors can be 
evaluated using the Cayleigh formula, as done again in [20] (formula B.l). We obtain 



i=i 


U(x)\ l l 1 - 1 

IP 


(3.21) 


which is indeed the Taylor expansion of the tree function. 

Actually, this direct derivation of the leading order provides us a nice interpretation of next-to-leading contribu¬ 
tion, the last of (3.17) 


(i) _ T(U{x)/2) (o) 

T i-T(U{x)/2y YT 


(x) = \u{x)T'(U(x)/2)VYp\x) = J2 


i=i 


‘ _vy (0) 

i\ T 


(x) 


(3.22) 


where, for the last equality, we have used the above Taylor expansion (of the derivative) of the tree function. In 
fact, following back in the reverse order the steps to the derivation of the formula (3.21), we can engineer a cluster 
expansion formula for Yf, 1 (x), 


yAx) 
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n U ^ d 2 t n e 2 p{ct H -(j )j )yY^\x), (3.23) 

p m cij >e E( Af,j) 


where we denoted A ~ B for A = B + 0(e 2 ). The extra symmetry factor l corresponds to the possibilities to graft 
a leaf to the tree Aj, l on any of its vertices, i.e. 


Y± 1] (a:) ~ U ( x ) ^ e. 
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h) / £ P(<t> 
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(3.24) 


In this formula, the extra vertex 4> is grafted to the tree at the vertex <f> a via a p-link. This vertex bear the potential 
Vl£ (0) instead of U(0). But our formula can also be obtained in another way: consider a tree ; and select a 
link that will later correspond to p(6 — 0 n ). Removing this link, one obtains two clusters, one which is rooted by x 
called Tf (in black in figure 1), and a subtree of the original tree denoted 77f with a new root 0 where the link we 
removed used to end (in red on figure 1). To compute the ^-correction to this specific y>-link, we can approximate 
at first order each of the two clusters and the generating function associated to T 2 is Y^ ( 0 ). The correction from 
the p-link in blue appears in the integral 


p{4>a - 0)£' O) (0)^: = ^T°\(j) c 


) + = a 


+ -£2 \(/> a ) + 0 ( e 2), 


(3.25) 









where 4> a is any vertex of Tf. The second term in the RHS is interpreted as the subleading £ 2 -correction brought by 
the blue p-link. It only remains to compute at first order the integral associated to the vertices of Tf which exactly 
corresponds to the formula (3.24). 


3.3 Contributions from articulation links 


The £ 2 -corrections provided by the articulation links are similar to those obtained in (3.24) for the trees. The only 
difference is that now the articulation link splits the cluster into two generic clusters instead of two trees. The 
correction k'^ rt (x) is thus obtained by grafting to a rooted cluster a leaf with the potential VT s [° 0 | t that generates the 
second cluster: 
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upon recalling (3.6) and U(x) = qQ(x)e x ( x \ Transferring the potential to the root, we find (by also using // = l/l 
and (3.6)) 






1 U{x) 

2 1 - U{x) 
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1 short 
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(3.27) 


The col lection to the free energy T^hort presents an extra factor (l/l and) 1/2 due to the possibility to exchange the 
two sub-clusters in absence of a root 



II II p(* - MvyJ^,(x). (3.28) 

iev(Af) <ij>eE(A f) 


By factorisation of the potential, we recover the integral Ij, which, eventually, leads us to the correction 


41m = 


1 r dx_ (o) 

4 J 2in short 




(3.29) 


4 Prepotential at subleading order 

Upon plugging the previous result for Z s h 0r t into the first of (2.11), we can give an explicit path integral expression 
for the Nekrasov partition function 6 

Z - (exp ^ j U 2 (qQ{x)e x ( x) ^j ^ + \J lo g (l - qQ(x)e x(x) ) Vlog (l - qQ(x)e x(x) 'j , 

(4.1) 

where the term involving the operator V (whose action on the logarithm is well-defined only perturbatively in q, at 
any finite order) yields the exact subleading correction. 

In the previous expression we may re-introduce the initial kernel k(x — y) in place of its inverse t(x — y) 
via a ghost field p(x) (Hubbard-Stratonovich transform) 7 * ; moreover, to better see the correction to the Nekrasov- 
Shatashvili action (which indeed contains two fields), we shall perform a change of variable ip(x) = —X(x), so 

6 In principle this procedure could be carried on at any order in £ 2 - 

7 Although convenient, the introduction of the extra field p(x) is not strictly necessary. It is indeed possible to perform directly a semi- 

classical treatment of the quantum field X. The equations of motion associated to the variation of X reads, at leading order, upon performing 






that 


Z ~ 


DpDp exp 


(^r j k ( x - y)p{x)p(y)dxdy + J p{x)tp(x)dx + ^ j U 2 ( qQ(x)e v( - x ^ ^ 

+ \ / log i 1 ~ vl °s i 1 * * * * * - q Q(x)e~^) . 

(4.4) 


In the r.h.s., the term of the second line is subleading in e 2 and will be neglected in the derivation of the equations of 
motion. The remaining term defines the action Sq [p, ip] that gives 8 
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= ~p(x) - log (l - qQ{x)e ~' p ' (a:) ) = 0. 


j dz 

6 2 S 0 5 2 S 0 

S 2 S 0 

6 2 S 0 

5p(x)5p(z) 5tp(z)5<p(y) 

5(p(x)5p(z' / 

5p(z)5ip(y)_ 


(4.5) 


These equations allows to determine the classical fields by solving a TBA-NLIE, 

p{x) = J k(x - y) log (l - qQ{y)e~^ y) ^j 2mp(x) = - log (l - qQ(x)e~‘ p< ' x) ^ . 

Once the classical fields are know, i.e. on-shell, we can compute (minus) the Hessian kernel defined by 

H(x,y) = - [ dz 

and find 

H{X ' y) = S{x - y) -h k(x ~ y)A( - y) ’ AW = 1 

The same result could be obtained from the Hessian kernel with only one field X (4.3) upon convolution with the 
kernel k(x — y), which is the contribution to the measure from the introduction of the field p according to (2.8), (2.9). 
Note that on-shell we also have A = e 2?7r/) — 1. Eventually, taking into account both classical action and one-loop 
corrections, we find for the prepotential T = e 2 log Z in terms of the classical solutions to (4.6) 

^ =\ j p(x)ip{x)dx + J Li 2 ( qQ(x)e ~^ x) ) ^ + y£ 2 j p(x)Vp{x)dx -y logdet H(x, y) + 0{el ), (4.9) 

where, for instance, the logarithm of the determinant of H(x, y) enjoys a Fredholm infinite series formula 


(4.6) 


(4.7) 


(4.8) 


log det H(x,y) 


I n*<*)^ n k - &+0 ~ ^i)- 


71—1 


71 — 1 


i =1 


1=1 


(4-10) 


a convolution with the original kernel k(x), as the following TBA-NLIE 

X(x) + J fc (x - y) log (l - qQ{y)e x(v) ) ^ = 0. 

Deriving once more, the Hessian kernel reads 

<j 2 S[X] _ 1 _ _ gQ(x)e x ^ 

5X(x)5X(y) (2i7r) 2 ^ 2in ^ 1 — qQ{x)e x ^X) 

Once convoluted with the kernel k(x — y) this gives the same as the Hessian kernel with two fields p and <p as computed below. 

8 At this order, it is not necessary to compute the variation of the term involving the operator V. 


( 4 . 2 ) 


( 4 . 3 ) 
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5 Study of 17(1) Af = 2 SYM 


The Nekrasov instanton partition function of U{ 1) SYM without matter is particularly simple, 


z U(i) = exp 



(5.1) 


However, the evaluation of the each term contributing to the prepotential exhibits a non-trivial dependence on the 
gauge coupling, which is remarkably canceled out order-by-order so to provide the simple result given above. This 
theory is thus a very good candidate to test our formula (4.9) for the subleading corrections in £2 to the prepotential. 
Here we perform this verification up to the order 0(A 4 ), i.e. four vertices. In particular, this is a good test of our 
conjecture about the non-contribution of irreducible p-clusters at the origin of our formula for Z s h ort [Y]. 

The potential associated to vertices for pure 17(1) SYM reads 


Q(x) 


1 

(x - a)(x — a + £1 + £ 2 ) 


(5.2) 


Note that the Coulomb branch vev a can be set to zero due to the invariance under translations of the instantons 
position. Furthermore, at the order of interest q and A can be identified since q = A + O (e 2 )- The field (f(x) is 
determined at the order 0(q 3 ) and 0(e 2 ) after solving by iterations the TBA-NLIE (4.6). Plugging the result in the 
second equation produces p(x) at the order ()(q l ). Isolating the singularities in the upper half plane at x = a + ne 1 
with n = 0,1, 2, 3 that contributes to the singular part p s j n „ , it is possible to compute Vp(x). From these expressions 
we deduce the following contributions after evaluating the integrals as sum over residues, 


1 f , w ^ q 2 , Q 3 4 , ( 119 2 49q 3 , 1583q 4 ^ , ^ 5 2 , 

2 J = - 2 ^f + I 24 - 3fcl + ( 84 - 72^ + Ssfef) ^ '**>’ 


Lh ( qQ{x)e 


~<p(x) 


+ 


dx q q 2 

2^ ~ + 2 ef _ 124 ' 364 

^ (q*_ _ 54 22 ^\ _ , „ (J> 


+-o — 


5 q 2 5 q 3 371q 4 \ 


4e 4 


+ 


84~m4^) E2 + 0W ' £2) ' 


l E ‘l J p(x)Vp(x)dx = - 4 + 27 1: ) + 0I '£4>- 


£2 , A \ ( q 5 q 2 49 q 3 2915 q 4 \ 5 2 ^ 

- - log det H(x,y) = ^ ^ ^ ^ j e 2 + 0 („ , e 2 ). 


(5.3) 


The determinant of the Hessian kernel is obtained by a consistent truncation of the Fredholm expansion (4.10) at the 
order n = 4. Taking the sum of these terms, we match (5.1) with J~u(\) = £2 log-Zjyn) = A/ei + 0( A 5 , e 2 )- 


6 Perspectives 

We have proposed a formula for the subleading corrections to the prepotential of SU(N c ) J\f = 2 SYM in the 
Nekrasov-Shatashvili limit. However, two physically justified assumptions are used for the derivation and shall 
deserve deeper analysis in the near future. First, the field theory argument rest upon the treatment of the exponential 
of the field as a meromorphic potential. This assumption can be bypassed by a treatment of the full partition function 
by Mayer expansion. Since this derivation is rather lengthy, we decided to present it elsewhere [24]. The second 
point that requires a closer examination, in the close future, is the e 2 -expansion of Z^ on from which the full partition 
function can now be derived. 

The results presented here are very general and may find application in numerous problems. For instance, the 
computation of planar amplitudes in J\f = 4 SYM involves series of nested integrals that are similar to those defining 
the Nekrasov instanton partition function [27]. In particular, they also exhibit at strong coupling a phenomenon 
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where some poles pinch the integration contours, which is responsible for the formation of a heavy ’meson’ and its 
bound states [28]. Hopefully, the formulae presented here might be generalised for computing the one-loop correction 
to classical minimal area result, namely the contribution at order of the inverse string tension £2 ~ 1 /V A (s/X is 
also the t’Hooft coupling). 

The presence of the TBA-like NLIE at leading order connects the supersymmetric gauge theories to the realm 
of quantum integrable systems. Turning on the £2 parameter defines an integrable deformation of these systems 
that is characterised by the presence of the SHc algebra. However, this deformation is still poorly understood in the 
perspective of integrable theories. With the derivation of the first order correction, we hope to clarify the situation and 
find a physical interpretation of the integrable deformation. It is noted that in the Nekrasov-Shatashvili dictionary £2 
plays the role of the inverse length of the thermodynamic limit. It suggests to interpret the £ 2 -corrections as further 
finite size corrections to the Yang-Yang functional. It would be interesting to find this interpretation in a simple 
integrable system, also for the development of the integrability theory. 

The possible connections with Painleve sixth equations and isomonodromy problems also deserve further inves¬ 
tigation [32, 33]. In the case of an SU (2) gauge group, the NS prepotential is related under the AGT correspondence 
to the semi-classical conformal block of Liouville field theory. This conformal block can be obtained by solving 
in the classical limit b —? 0 the null vector decoupling equation (NVDE) obeyed by a five-points correlator with a 
degenerate operator of level two inserted. This operation corresponds to insert a surface operator in the AGT dual 
instanton partition function [34], The NVDE is identified in [32] with a Schrodinger equation obtained in a quan¬ 
tised version of the sixth Painleve equation [35]. Our result should be connected to the 0(K) corrections to the WKB 
solution of this Schrodinger problem. 

Still in the case of SU (2) SYM, it is possible to take the colliding limit of the potential Q(x) in which the leading 
order prepotential compute the classical irregular conformal blocks [36]. This quantity is related to the accessory 
parameter of a Mathieu equation (i.e. the energy in the Schrodinger form) obtained as the aforementioned classical 
limit of the NVDE [37], The formula presented here allows the computation of the first semi-classical correction to 
the irregular conformal block. This correction could be related to a deformation of the Mathieu equation. 

Acknowledgements. It is a pleasure to thank A. Lerda and R. Tateo for useful discussions. This project was 
partially supported by INFN grant GAST, the UniTo-SanPaolo research grant Nr TO-Call3-2012-0088, the ESF 
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